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The processes leading to nonthermal condensate vaporization and charge-density wave (CDW) 
melting with femtosecond laser pulses is systematically investigated in different materials. We find 
that vaporization is relatively slow {jv ~ 1 ps) and inefficient in superconductors, exhibiting a strong 
systematic dependence of the vaporization energy J7„ on T^. In contrast, melting of CDW order 
proceeds rapidly (xm = 50 ~ 200 fs) and more efficiently. A quantitative model describing the 
observed systematic behavior in superconductors is proposed based on a phonon-mediated quasi- 
particle (QP) bottleneck mechanism. In contrast, Fermi surface disruption by hot QPs is proposed 
to be responsible for CDW state melting. 



Photoinduced phase transitions on the sub-picosecond 
timescale in superconductors and other electronically or- 
dered systems have attracted increasing attention in re- 
cent years, partly because of the fundamental desire to 
control collective states of matter with femtosecond laser 
pulses, partly because of potential applications in ul- 
trafast phase-change memories, and partly because they 
may reveal some fundamental insight into the mechanism 
of high-temperature superconductivity. However, the de- 
tails of how absorbed photons cause a change of state 
on ultrafast timescales have so far not been investigated 
in detail. Recently, a study on La2_£cSr^Cu04 (LSCO) 
single crystals by Kusar et al^ showed that the energy 
required to vaporise the superconducting (SC) conden- 
sate Uy can be determined with reasonable accuracy. 
These and other measurements on cuprate superconduc- 
tors since thcEP® give values of C/„ which are sometimes 
significantly larger than the experimental condensation 
energy C/c, or the BCS theoretical values^, opening the 
question of mechanism for the destruction of the con- 
densate, the dependence of on doping and on the 
critical temperature Tc. It is also not clear if this oc- 
casional large discrepancy between Uc and is peculiar 
feature of cuprate superconductors, is a consequence of 
a large gap in the density of states, or depends on the 
detailed mechanism responsible for the formation of the 
low-temperature ordered state. 

To try and answer these questions, we present the first 
systematic study of the dependence of on doping and 
Tc for cuprates and compare our data with the iron pnic- 
tides (Ba(Feo.93Coo.o7)2As2 and SmFeAsOo.8Fo.2), the 
conventional superconductor NbN and two large-gapped 
charge-density wave systems K0.3M0O3 and TbTca. We 
made a special effort to determine the geometrical factors 
and optical constants as accurately as possible in order 
to achieve the best possible accuracy in determining Uy. 
We find that the observed systematic behavior reveals 
the mechanism for breaking up the collective state in 
high- Tc superconductors to be relatively inefficient and 
phonon dominated. This is in stark contrast to CDW 
systems, where the electronic destruction mechanism is 



faster, more direct and efficient. We present a model 
to describe the observed systematics in the cuprates and 
discuss the markedly different mechanisms for the pho- 
tonic destruction of the low-temperature ordered state in 
the two types of system. 

Our experiments were performed using a standard 
pump-probe technique with a 800 nm laser with 50 fs 
pulses as described in ref.^^. A probe wavelength of 
800 nm is chosen because it is a well-understood probe 
wavelength for investigating the superconducting (SC) 
and pseudogap (PC) response in pump-probe experi- 
ments and the best signal/noise ratio as determined 
by previous studies. The 2-component response, PC 
and SC, including lifetimes, T-dependence, anisotropy 
and doping dependence have been previously studied at 
800 nm and other pump/probe wavelengths as well, in- 
cluding pump/probing in the THz region. The results 
of these measurements are quantitatively self-consistent, 
and have all been systematically discussed. The second 
reason for using 800 nm (1.5 eV) is that the optical con- 
stants are known more accurately at this wavelength bet- 
ter than any other, which improves the accuracy of the 
determination of the vaporization energy. The reason for 
not using far infrared (FIR) wavelengths in the gap re- 
gion is that even with transform-limited pulses, the pulses 
become too long to be able to distinctly detect the SC 
component and distinguish it from the PC. Hence THz or 
FIR pulses in the gap region are not suitable for our pur- 
pose. Recently the electron distribution function ther- 
malization has been measured to be ^ 60 fs in LaSCO 
and - 100 fs in YBCCPil, while the quasiparticle (QP) 
signal used for detection of the presence of the supercon- 
ducting state in the cuprates has a characteristic lifetime 
nearly two orders of magnitude longer and can thus be 
easily distinguished from the PC and energy relaxation 
which are both much faster. In order to accurately deter- 
mine the deposited energy, pump and probe laser beam 
diameters were accurately measured with calibrated pin- 
holes. The light penetration depth Aop and reflectivity 
R, are obtained from published optical conductivity and 
dielectric function data in each casd^. The measurements 
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were performed on freshly cleaved YBa2Cu307_5 single 
crystals with TcS of 90, 63 and 60 K, respectively and 
Ca^Yi_^Ba2Cu307-5 with x = 0.132 and 5 = 0.072 
and Tc = 73 K. The photoinduced reflectivity transient 
Ai?/i? for YBaaCugOe.s well below (4 K) and above 
Tc (at 68 K) are shown in Fig. la), while the net super- 
conducting signal, obtained by subtracting the response 
at T = 68 K from the response at T = 4 K : AR.JR 
= AR/ R\t<t,-AR/ R\t>t, is plotted in Fig. lb) for 
different laser pump pulse fluences. This behavior is 
generic, so we have shown the raw data only for one sam- 
ple. ARsc/R shows a saturation plateau above > 100 
/iJ/cm^, which is a signature of the destruction of the 
SC state. 

The normalized magnitude of ARsc/ R for different ma- 
terials and different doping levels as a function of is 
shown in Fig. 2. At low fluence ARgc/R is linearly 
dependent on up to a certain threshold fluence J-t, 
whereupon it either becomes saturated or continues in- 
creasing, but with a different slope. To determine the 
vaporization threshold fluence from the data in Fig. 
2, we use the inhomogeneous excitation model^. To take 
into account carrier-density dependent change of AR/R 
arising from e — h asymmetry in the band structure, we 
have also included a linear term in the fit to the data in 
Fig. 2, such that AR/AR^ = ARsc/AR, + Em- 
pirically, we find that in LSCO, </> ~ while in YBCO 
samples ~ 4>q^, where (fiQ is a constant. Note that 

this term has no effect on the values of Ft obtained from 
the fits to the data in Fig. 2. (As a check, Fj- can also 
be directly read out as the point where ^^'^(■^'> departs 
from linearity). 

A plot of the vaporisation energy [1^ ~ {I — R)Ft/Xop 
obtained for cuprates, the iron pnictides and NbN is 
shown in Fig. 3. We observe systematic increase of [/„ 
with Tc, giving an approximate square power law depen- 
dence Uv = rjTci with 77 ~ 1.5 x lO^^K^^atom^^ irre- 
spective of whether the material is overdoped, optimally 
doped or underdoped. 

In Table 1 we compare with the condensation en- 
ergy [/"''P obtained from specific heat measurements. In 
the cuprates, we notice a significant discrepancy between 
the magnitudes of [/„ and f/^^^. While U^^^ shows a 
distinct asymmetry between overdoped and underdoped 
materials with the same T^^, Uy depends systematically 
on Tc irrespective of doping level. Thus [/„ and U^^'^ 
show fundamentally different systematics. 

To understand this apparent discrepancy between the 
condensation and vaporisation energy, let us now con- 
sider the mechanism of energy transfer between photoex- 
cited (PE) carriers and the condensate. After absorption 
of a photon, the PE electrons and holes lose energy very 
rapidly through "avalanche" scattering with phonons on 
a timescale of te < 100 fs for YBCO and < 60 fs for 
creating a large non-equilibrium phonon popu- 
lation in the process. The phonons whose energy exceeds 
the gap hujpi^ > 2A can subsequently excite QPs from 
the condensate, but these recombine again and so a bot- 
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Figure 1. (Color online), a) AR/R for YBaaCugOe.s for 
T<Tc (squares) and for T>Tc (circles) at excitation fluence 
F = 17.4 /xJ/cm^. b) ARsc/R for YBa2Cu306.5 as a func- 
tion of time delay for different F, showing saturation above 
46/1 J/cm^. The arrow indicates the direction of increasing 
F. The condensate vaporization time r„ is between 0.5 and 
1 ps. (At high fluences the PG signal interferes with the 
measurement, so t„ cannot be more accurately measured in 
YBaCuO). 

tleneck occurs, in which high-frequency phonons are in 
temporarily in quasi-equilibrium with QPs as described 
by Kusar et al. for LaSrCuC^. The process eventually 
causes the destruction of the condensate within 0.5 ^ 1 
pfP. 

The phonons with huj > 2A can contribute to 
the breakup of the condensate, but the low frequency 
phonons with Hui < 2A also created in the initial 
avalanche do not have enough energy to excite QPs, so 
do not contribute to the vaporisation process. This lost 
energy to low-energy phonons is: 

[/lost = / Af{uj)D{uj)hojdoj (1) 
Jo 

where D{ijj) is the phonon density of states and A/(w) = 
/Ar£;(w) — fsi^) is the difference between the non- 
equilibrium and equilibrium phonon distribution func- 
tions. A/(a;) is material-dependent and depends on the 
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Material 


Ta (K) 


J-'t (m J/cni^ ) 


Aop(mn) 








;7v/(7°^P 


T„ (ps) 


Lai.85Sro.i5Cu04 (opt)^'- 


38.5 


5.8±2.3 


150 


2.6±1.0 K/Cu 


0.3 K/CJD 




> 8.5 


0.9 


YBaaCuaOe.s (ud) 


60 


13.1±0.5 


90 


8.3±Q.3 K/Cu 


0.38 K/CiJSl 




22 


0.5~ 1 


YBaaCuaOe 6 (ud) 


63 


12.1±1.3 


85 


8.1±Q.9 K/Cu 


0.45 K/CiPl 




18 




YBaaCuaOr (opt) 


90 


10.8±0.6 


66 


9.3 ±0.5 K/Cu 


1.9 K/CiPl 
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Yo.sCao,2Ba2Cu306.9 (od) 


73 


14.3±1 


75 


10.8±0.8 K/Cu 


1.4K/CJS1 




7.7 




BiaSrsCaCuaOs+i (udP 


78 


16±1 


124 


9±0.8 K/Cu 


1.8 
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SmFeAsOo.sFo.J^^ 


49.5 


2.6±0.2 


31 


1.8±0.1 K/Fo 


- 




- 


0.3-0.5 


Ba(Feo.93Coo.o7)2As2 


23 


0.43±0.04 


34 


0.3±0.03 K/Fc 


0.15±0.02' 




2 


0.5 


NbrJS 


16 


t/v = 23 mj/cm^ 




0.24 K/Nb 


0.14K/Nb 




1.7 






Ta (K) 


J-'t (M'Vcni^ ) 








jjthcory 


C/v/(7„ 




TbTe;^i^ 


336 


47±2 




52±2 K/Tb 




4Q.6±1.7 K/Tb 


1.3 


o.iia 


Ko.aMoOipHIl 


180 


105±17.5 




35±10K/Mo 


>10K/Mo 


30±10 K/MgW 


1.1 ~ 3.5 


0.05-0. fiS 



Table I. Tc, J^t, Uv and Uc for different materials. Uc and (7v are calculated per metal atom of Cu (planar), Fe, Nb, Tb or Mo. 
For LSCO and YBCO, the data for Uc are maximum experimentally measured value^^. S in YBCO was changed by annealing 
in flowing oxygen, while TcS were determined from magnetization measurements. ForTbTes and K0.3M0O3, the electronic 
part of Ue of the condensation energy is estimated using available values of A and 

jYjmH!. For YBa2Cu30.^_5, the penetration 
depths arfPl 66±15 nm {5 = 0), 85±15 nm {5 = 0.4), 90±15 nm {S = 0.5), and 75±15 nm for CaYBCO. 
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Figure 2. (Color online). A-Rsc as a function of normalised 
fluence T (Tt- Note the linear response below Tt ■ 

electron phonon coupling, more precisely on the Eliash- 
berg coupling function o?F{lS), but is not a very strong 
function of cJ^^l, Xo estimate the dependence of [/;ost 
on Tc from Eq. (1), we assume that A/(a;) is constant, 
and approximate the experimental D(u}) by an inverted 
parabola D{uj) ~ ^uj{2ujo — oj), where a is a constant. 
Integrating Eq. (1), we obtain: 

f/iost^aA^l'^-f]. (2) 

Assuming a constant gap ratio 2A/A:b7c = R, with R = 
4, for YBCO, hujQ = 40 meVl22l extending D{uj) to 80 meV 
as shown in the insert to Fig. 3. Eq. (2) gives the curve 
shown in Fig. 3. For iron pnictides, Hujq ~ 2(j221, but 
the predicted variation of Uiost on Tc is not significantly 
different. We see that for the superconductor series Eq. 
2 predicts the dependence of Uy on Tc very well, which 
is not surprising, considering the gross features of D{uj) 
do not vary significantly. 




Figure 3. (Color online). AUv expressed in K per planar 
Cu as a function of Tc^ for the cuprates. The data for NbN, 
SmFeAsOo.8Fo.2 and Ba(Feo.93,Coo.o7)2As2 are included for 
comparison (in K/Fe or K/Nb). The solid curve is a plot of 
l7iost using D{ijjph) appropriate from Eq. (2) for YBaCuO 
(red) and the iron oxy-pnictide (blue). The dashed line is a 
square law Uv ~ ■ The insert shows the phonons with 
huj > 2A (shaded) which can break pairs. The measured 
phonon density of states -D(tjph) for YBCC^'' is approximated 
by a parabola (line). 



The total vaporisation energy is the sum of the conden- 
sation energy and the lost energy, Uv — Uc + U\ost- Since 
for large gap systems Uy » 11°^^ we have Uy ~ ?7iost, so 
U^^P is thus only a small contribution to Uy explaining 
why the anomalous doping dependence of [/^^p in the 
cuprates is not displayed by Uy. 

Comparing Uy for NbN measured using the same tech- 
nique (Table I), U^/U^^^ = 1.7 is considerably smaller 
than in the cuprates, which can now be understood in 
terms of Eq. (1) and the insert to Fig. 3: When 
2 A << ^WDebye almost all phonons can excite QPs, so 
Uiost 0, and Uy ~ Uc, in agreement with the data. For 
small-gap superconductors, the optical method can be 
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used to estimate the superconducting condensation en- 
ergy. The data on the two iron pnictides Hsted in Table 
I also seem to follow the predicted behaviour in Fig. 3. 

We now turn our attention to the data obtained us- 
ing the same technique on two different charge-density- 
wave materials TbTes and K0.3M0O3. Both have 
CDW gaps which are much larger than any phonon 
energy: 2/!^cd w ^ f^D (Act)IV=125 meV and 75 
meV respectiveljCmil) . K0.3M0O3 is considered a pro- 
totypical Peierls system with a full gap in the elec- 
tronic density of states opening at the metal-insulator 
transition^^' Moreover, in KM00.3O3 the electron-phonon 
coupling constant {\e-p — 0.35) is comparable with the 
cuprates {\e-p — 0.2 5 for YBaCuO and Ag-p — 0.5 
for LaSCO respectively ^^P^III] go with a phonon-mediated 
CDW melting mechanism we would expect a very large 
Um/Uc ratio. Instead, it is significantly smaller than 
in the large-gapped cuprates, with Um/U^^^ ~ 3.5. In 
TbTe3 the ordered state opens a partial gap only along 
certain directions on the Fermi surface, the rest remain- 
ing ungappecP^. Thus, if the same phonon-mediated 
vaporisation mechanism were operative, the ungapped 
CDW system should behave very differently to KM00.3O3 
on one hand, and the superconductors on the other. 
Comparing Uv with the electr onic contribution to the 
CDW condensation energjESES = Nq/S.'^ {\+log'^) , 
where is the Fermi energjEIl we see in Table I that the 
predicted ratios are Um/Uc ~ 1.3 for TbTes and 1.1 for 
K0.3M0O3, confirming that Uiost and implying that 
phonons are not significantly involved in the CDW gap 
destruction and CDW melting process. Thus, in spite 
of similar e-p interaction strength and very large gaps, 
Um/Uc is consistently smaller in the CDWs than U^/Uc 
in cuprates, highlighting the very different mechanism for 
the destruction of the low-temperature CDW state. 

The reason is fundametal: Unlike in superconduc- 
tors, PE carriers in CDW systems can cause a direct 
rapid disturbance of the Fermi surfac^l^ on a timescale 
shorter than the QP recombination time. It is known 



experimentally'^ and theoreticallj^that the CDW state 
is extremely sensitive to charge imbalance, which dis- 
rupts the Fermi surface (FS) nesting. In photoexcitation 
experiments, overall charge neutrality is preserved, but 
immediately after photon absorption, even the slightest 
electron-hole band asymmetry causes a transient shift of 
the chemical potential, causing a disturbance 6g of the 
FS. This causes a destruction of the nesting at kp and a 
suppression of the divergence of the electronic suscepti- 
bility at 2kF c ausin g the electronic charge density wave 
to rapidly meliP^fil. A fast destruction time in CDWs is 
experimentally confirmed: in K0.3M0O3, t„i = 50 ^ 100 
fgiSl, which is an order of magnitude less than Ty ^ 0.9 ps 
in LaSCO, which has a similar e-p interaction strength. 
Indeed, quite generally, the melying time is signifi- 
cantly shorter than the in superconductors, as shown 
in Table I (not all risetimes are reliably measured so 
far, so only the unambiguous ones are listed). With 
T„i = 50 ^ 100 fs, the ions can just barely move into new 
equilibrium positions within this time, (the l/4^eriod of 
the amplitude mode in K0.3M0O3 is ~ 150 fs)W 

We conclude that photon-induced change of state in 
superconductors and CDW systems relies on very dif- 
ferent mechanisms, occuring on different timescales and 
requiring different amounts of energy. Systematic and 
experimentally precise measurements of the dependence 
of Uy on doping and on Tc in cuprates, pnictides, low- 
Tc superconductors show that vaporisation can be de- 
scribed by the established Rothwarf- Taylor QP recombi- 
nation mechanism, in a process where energy transfer to 
the condensate is slowed down and made inefficient by 
the phonon-QP relaxation bottleneck. In CDW systems 
on the other hand, the dominant mechanism for the de- 
struction of the ordered state is electronic, so in spite of 
the fact that the CDW gaps 2Acdw ^re typically much 
larger than the highest phonon frequencies, and electron- 
phonon coupling strengths are similar as in cuprates, the 
destruction is faster and more direct. 
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